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Abstract. 

The interstellar medium (ISM) in galaxies is multiphase and cloudy, with stars forming in 
the very dense, cold gas found in Giant Molecular Clouds (GMCs). Simulating the evolution of 
an entire galaxy, however, is a computational problem which covers many orders of magnitude, 
so many simulations cannot reach densities high enough or temperatures low enough to resolve 
this multiphase nature. Therefore, the formation of GMCs is not captured and the resulting 
gas distribution is smooth, contrary to observations. We investigate how star formation (SF) 

0~ • proceeds in simulated galaxies when we obtain parsec-scale resolution and more successfully 
capture the multiphase ISM. Both major mergers and the accretion of cold gas via filaments are 
dominant contributors to a galaxy's total stellar budget and we examine SF at high resolution 
y—n 1 in both of these contexts. 
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By performing adaptive mesh refinement (AMR) simulations in which we resolve the 
multiphase ISM and the formation of GMCs, we investigate the nature of galactic SF in 
two important contexts: major mergers (Section 1) and cold accretion via filaments at 
_-, . high-redshift (Section 2). 
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OO I 1- Star formation modes during major mergers 

To date, we have performed 8 simulations of 1:1 galaxy mergers (M ga j « 8 x 1O 1O M0) 
with a maximum 5pc spatial resolution and 1.6 x 10 4 M Q effective gas mass resolution 
using the AMR code RAMSES (|Tevssierll2002t ) , providing us with a representative sample 
that contains a variety of orbital parameters. To achieve our high resolution, we use a 



o 

'pseudo-coo ling' equation of sta te which dictates the temperature of gas based on its 
density (see iTevssier et aL 20101 for details). This approach has been shown to produce 



an ISM density power spectrum that agrees well with observations (Bournaud, 2010 in 
•rH . prep.). We choose a moderately high density threshold (~ 10 3 atoms cm -3 ) and a very 
low efficiency (~ 10~ 4 ) for t he SF prescription, and e mploy the Sedov blastwave model 
of supernova (SN) feedback ()Dubois fc Tevssier1l2008f) . This choice of parameters gives 
a star formation rate (SFR) of ~ lM Q yr _1 in the galaxies when they are simulated in 
isolation, which is reasonable for az = disc galaxy. Work is currently underway to test 
the sensitivity of the SFR to changes in these parameters. 

1.1. Extended star formation 

In the classical picture, merger-induced starbursts result when tidal torques drive large 
amounts of gas into the centre of the merging structure, where it is compressed and forms 

1 



2 



Leila C. Powell et al 




Figure 1. Left and Middle: Gas density images of the merging galaxies, one face-on, one 
edge-on (left) and of one of the same galaxies when evolved in isolation (middle) at t=1.5 
Gyr. Red contours indicate the regions with densities above the SF threshold. Size of image is 
~ 13kpc. Right: SFR versus time for both of the merging galaxies (solid line) and the sum of 
the SFRs in the same two galaxies when evolved in isolation (dashed line). 



stars rapidly (Struck, this meeting). We report the occurrence of another mode of merger- 
induced SF: clumpy, extended SF. In Fig. Q] we clearly demonstrate this alternative 
mode by comparing the gas density map of the merging galaxies (left) with that of 
one of the same galaxies when evolved in isolation (middle) at the same time instant 
(t = 1.5 Gyr). Red contours indicate the high density, star- forming clumps and these are 
considerably more numerous and significantly less centrally concentrated in the merging 
galaxies (where the clumps are spread over a region of ~ 7kpc in diameter) than in the 
isolated galaxy (clumps spread over « 3kpc). This is particularly clear when comparing 
the merging galaxy seen face-on in the left-hand panel of Fig. [1] with the isolated galaxy 
(middle panel). By reference to the time evolution of the SFR (Fig. Q] right), we can see 
that t — 1.5 Gyr occurs just before the maximum SFR, indicating that this extended, 
clumpy SF is responsible for the starburst which has increased the SFR in the merger 
(solid line) by a factor of ~ 10 relative to the SFR in the same two galaxies evolved in 
isolation (dashed line) i.e. the starburst occurs before the final coalescence of the two 
galaxy centres. 

This type of merger-induced st ar cluster formation over an extended region is also 
observed in the Antennae galaxies. Tevssier et al. (|2010() have demonstrated that this is 
only reproduced in a simulation of the Antennae with RAMSES if a resolution of ~ lOpc 
is achieved; they find that the star formation is no longer clustered if the gas is thermally 
supported (as occurs at lower resolutions), highlighting the strong dependence of the 
mode of simulated SF on resolution. 

In Fig. [5] we show some examples of the other varied star-forming features found 
scrcndipitously in our merger sample. 




Figure 2. Gas den- 
sity images from two 
different mergers. 
Left: Formation of 
'beads on a string' 
of length ~ 100 kpc. 
Right: Two tidal 
dwarf galaxies (one 
with M ~ 10 9 M Q ) 
formed after the first 
pericenter passage. 
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2. Star formation during cold-mode accretion and the onset of a 
galactic wind 

The Nut Simulation Q is an ultra-high resolution cosmological resimulation o f a Milky 
Way l ike galaxy forming at th e intersection of 3 filaments in a ACDM cosm ology (ISlyz et al 
l2009t iDevriendt et al.lEoioj ). We per form the simulation wi th RAMSES (jTevssierl l2002h 



including cooling, a UV background (Haardt fc Madaulll996l ). SF and SN feedback with 



20081 ) achieving a maximum physical spatial reso- 



metal enrichment ([Dubois fc Tevssierl ll 
lution of w 0.5pc in the densest regions at all times. A thin, rapidly rotating disc forms, 
which is continuously fed by cold gas that streams along the filaments at ~Mach 5. The 
disc is gravitationally unstable and fragments into star-forming clumps (see Devriendt 
et al (2010, in prep.) for details of the disc and Slyz et al (2010, in prep.) for a study 
of the clumpy star formation). SF in the galaxy, its satellites and haloes embedded in 
the filaments results in many SN explosions, which are individually resolved. The numer- 
ous bubbles overlap, creating an extended ca vity through which the hot gas escapes the 
galactic potential in the form of a wind (see ISlyz et al 1 120091: IDevriendt et al.|[2010l for 
figures illustrating this and additional technical details). By z = 9, when the dark matter 
halo mass has reached ~ 10 9 Mq, the wind extends to ~ 6r V i r , filling the virial sphere 
with hot, diffuse gas not unlike the shock-heate d gaseous halo found in galaxies above 
the mass threshold for a stable virial shock (e.g. iBirnboim fc Dekelll2003r ). To date, the 
simulation has reached z « 7.5 and it is still running. 

Here we present results from a study of the interaction between the cold filaments and 
the hot, far-reaching galactic wind and examine the impact this competition has on the 
SFR (Powell et al, 2010 in prep.). 



2.1. Inflows and outflows: measuring accretion and winds 

To assess the relative importance of the filaments and the galactic wind we measure the 
mass inflow and outflow rates within the virial radius of the galaxy. We divide the gas into 
various phases with the use of temperature and density threshol ds, similar to the ap proach 
in other simulation studies of hot and cold mode accretion fe.g. lAgertz et al.ll2009D . Note, 
however, that the thresholds are resolution-dependent so ours differ from those in the 
previous lower resolution studies. The categories we define are: hot diffuse (T > 2 x 10 5 K, 
n < 10 atoms cm -3 ), warm diffuse (2 x 10 4 K < T < 2 x 10 5 K, n < 10 atoms cm -3 ), cold 
diffuse (T < 2xl0 4 K, n < 0.1 atoms cm -3 ), filaments (T < 2xl0 4 K, 0.1 < n < 10 atoms 
cm -3 ) and 'clumpy' (n > 10 atoms cm -3 ). We verify that these categories accurately 
separate the different gas components using 3D visualisation software. 

In Fig. [3]we compare the inflow in a control run with no SN feedback (left) with that 
in the Nut simulation (middle) and examine the outflow in the Nut simulation (right). 
We highlight several key results: 1) the inflow of cold gas in both runs (first two panels), 
is dominated at all radii by the contribution from the filaments (blue dotted line) , which 
supply material right down to the disc, 2) the total inflow rate (first two panels, black solid 
line) is almost identical in the simulations with SN feedback (middle, the Nut simulation) 
and without SN feedback (left, control run), revealing that the galactic wind has had a 
negligible effect on the accretion rate, 3) the only significant mass outflow occurs in the 
Nut simulation (right) yet this is only around 10 per cent of the mass inflow rate and 4) 
there is a non-negligible cold diffuse component (blue dashed line) in this outflow, which 
is gas that was previously inbetween the filaments that has been entrained by the wind. 



t Nut is the Egyptian goddess of the night sky 
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Figure 3. Gas inflow averaged in 100 pc physical spherical shells out to r v ir in the control run 
without supernovae feedback (left) and the Nut simulation (middle) at z = 9. Gas outflow in 
the Nut simulation (right). Flow rate of all gas (solid black line), subdivided into dense (blue 
dot-dash line), filamentary (blue dotted line), cold diffuse (blue dashed line), hot diffuse (red 
dashed line), warm diffuse (orange dashed line) and hot+warm diffuse (red solid line). 

2.2. What happens to the star formation rate? 

We have established that the high-velocity, far-reaching galactic wind has been unable 
to impact the supply of cold gas, the 'fuel' for SF, because it is delivered to the galaxy 
in supersonic, highly collimated, dense filaments. The other major influence on the SFR 
is the efficiency with which this fuel is converted into stars. 
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Figure 4. SFR of main 
galaxy averaged over lOMyr 
Left: SFR divided by cold 
gas inflow rate for control 
run (blue) and Nut simulation 
(red). Right: SFR in Nut 
simulation divided by SFR in 
control run. Horizontal lines 
indicate when ratios = 1. 



The left-hand panel of Fig. [4] shows the ratio of the SFR over the net cold gas inflow 
rate for the control run (blue) and the Nut simulation (red). For both runs the ratio 
oscillates around 1 (the horizontal line), indicating that the conversion of the filament- 
supplied cold gas to stars is very efficient. The mean efficiency in the central star- forming 
region is ~ 0.1 in both cases; an order of magnitude higher than the 1 per cent efficiency 
used in RAMSES to dictate the SFR on a cell by cell basis. The right-hand panel of 
Fig. H] shows the ratio of the SFR in the Nut Simulation to that in the control run (the 
actual SFR is a few Moyr -1 ). Approaching z = 9, the SFR in the Nut simulation (with 
SN feedback) exceeds that in the control run (without SN feedback). We attribute this 
positive feedback effect to the metal enrichment in the Nut simulation (absent in the 
control run), which facilitates cooling of the gas, boosting SF. 
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